Bow wave from ultra-intense electromagnetic pulses in plasmas 
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We show a new effect of the bow wave excitation by an intense short laser pulse propagating in 
underdense plasma. Due to spreading of the laser pulse energy in transverse direction, the bow wave 
causes a large-scale transverse modulation of the electron density. This can significantly increase 
the electric potential of the wake wave since the wake wave is generated in the region much wider 
than the laser pulse waist. 

PACS numbers: 52.38.-r, 52.38.Kd, 52.65.Rr 

Bow waves are observed at various conditions in fluids, gases and plasmas. A bow wave is formed in the water at 
the bow of a moving ship [1^ in addition to Kelvin waves left behind; it defines the outer boundary of the ship wakes. 
Reducing this effect is an important aim of shipbuilding, since bow waves take away a portion of the ship energy 
and can damage shore facilities. In air, a spectacular bow wave accompanied with a vapour cone manifesting the 
Prandtl-Glauert singularity [2\ is excited by an aircraft when it breaks the sound barrier. In space, plasma bow waves 
appear at various scales. A bow shock is formed where the Earth's magnetosphere is squeezed by the solar wind flow 
[3] . The largest ever seen bow waves were observed in the collision of galaxies by the Spitzer Space Telescope [4] . As 
these examples show, bow waves lead to a vast range of physical processes such as a transverse transfer of momentum 
and matter and particle acceleration. 

As is well known, an intense short laser pulse propagating in underdense plasma excites wake waves, similarly to a 
moving ship producing the Kelvin waves. This analogy underlying a fruitful exchange of ideas between hydrodynamics 
and plasma physics can be extended further. In this Letter we demonstrate that a laser pulse can excite a bow wave 
in coUisionless plasma. Like a bow wave from a ship, the bow wave from a laser pulse spread aside the energy and 
the momentum of the laser pulse. Although generation of shock waves in laser plasmas has been known shortly after 
the flrst laser was built [5], to our best knowledge, the bow waves formed by relativistic electrons at the ultra-short 
laser pulse propagating in plasma have not been discussed so far. 

The laser-driven wake waves in plasmas are used in the laser wake field accelerator (LWEA) [6^ and the flying mirror 
(FM) [7j concepts. A tightly focused and sufficiently intense laser pulse excites a wake field |B] in the so-called blowout 
regime [H [9] or bubble regime [10], where, in addition to a longitudinal push, the laser pulse expells electrons also 
in transverse direction, forming a cavity void of electrons in the first period of the wake wave. Thus a finite waist of 
the laser pulse results in a transverse motion of electrons, which leads to a transverse wave-breaking |llj underlying 
one of the mechanisms of self-injection of electrons into the accelerating phase of the wake wave |S]. We show that 
when the waist of a sufficiently intense laser pulse is less than the wake wave length, the expelled electrons can not be 
confined near the cavity and their transverse outflow forms a bow wave. The bow wave should be distinctly seen in 
the limit of high intensity and narrow waist laser pulses and low-density plasma, which correspond to the conditions 
for generating higher energy electrons by laser wakefield acceleraion [8j [12] . 



*Also at A. M. Prokhorov Institute of General Physics of RAS, Moscow, Russia. 
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FIG. 1: (color). Electron density (coded into color, C, and opacity, O, with ray-tracing) in the wake of the laser pulse, which 
is represented by the electromagnetic energy density isosurface [E^ + B^^im^Locj ef' = 4 (red), at t = 100 x 2-k/ui. 

We present the results of three-dimensional (3D) particle-in-cell (PIC) simulations with the Relativistic Electro- 
Magnetic Particle-mesh (REMP) code [15 performed on Altix 3700 supercomputer at JAEA-Tokai. The laser 
pulse with the initial intensity / = 6 x lO^^W/cm^ x (1/im/A)^, corresponding to the dimensionless amplitude 
a = eEQ/umeC — 6.62, propagates along the x-axis; it is linearly polarized in the direction of the y-axis and has the 
size of £x X £y X £z = lOA x lOA x lOA. The pulse is incident on a fully ionized plasma slab with the electron density 
He — 1.14 X lO^^cm^"^ X (1/im/A)^. Here A and w are the laser wavelength and frequency, e and nie are the charge 
and mass of electron, Eq is the magnitude of the laser pulse electric field, and c is the speed of light in vacuum. 
We consider the electron bow wave, when the ion response can be neglected, i. e. the ion-to-electron mass ratio is 
mi/nie — > 00. The simulation grid dimensions are 4000 x 992 x 992 along x, y and z axes; the grid mesh sizes are 
da; = A/32, dy = dz = A/8; total number of quasi-particles is 2.3 x 10^°. For observing the bow wave, the number 
of quasi-particles must be sufficiently high since the density of quasi-particles expelled in transverse direction drops 
with a distance r from the axis as (in 3D configuration). In the figures, the time and spatial units are the laser 
wave period and wavelength, respectively. 

The laser pulse excites a strong wake wave, whose first period forms a cavity completely void of electrons. Figs, [l] 
[2] In addition to a longitudinal push, the laser pulse drives electrons also in transverse direction, forming the bow 
wave. As seen in Fig.|2ja), the longitudinal electric field (Ex) is formed in a region which is much greater in transverse 
direction than the cavity width corresponding to the first period of the wake wave. Due to this, the electric potential, 
(fi, is not limited by the transverse size of the cavity. We note that at present conditions the length of the cavity (i. e., 
its size along x-axis) is 1.5 times greater than its width (the size in transverse direction). In the regimes not favourable 
for the bow wave formation, e. g. when the bow wave does not detach from the cavity, the electric potential in the 
cavity is determined by the smallest dimension of the cavity, either longitudinal or transverse. 

A strong transverse wave-breaking (TWB) [11 forms the transverse outflow of electrons, which increases the electric 
potential in the subsequent periods of the wake wave, similarly to the bow wave. However, in contrast to the bow 
wave, which arises because electrons are directly pushed by ponderomotive force of the laser pulse, the TWB wave is 
due to an electrostatic field of the wake wave itself caused by the finite waist of the laser pulse. 

The electron motion in the wake wave, bow wave as well as motion due to TWB, forms a transient electron density 
modulations with associated electromagnetic field. Fig. [3] shows the magnetic field co- moving with the cavities. The 
polarity of the magnetic field corresponds to the forward motion of a negatively charged object. The strong magnetic 
field associated with the region of the electron density compression due to TWB can be described in the approximation 
of a Lienard-Wiechert potential of a moving point charge PJ] . 

As shown in the electron phase plot, Fig.|4]jc), the transverse momentum of the electrons undergone TWB behind 
the first cavity is higher than that of electrons in the bow wave. This is due to Coulomb explosion of the electron 
flow concentrated near the axis. Fig. |4][a,b). 

The shape of the cavity in the electron density behind the laser pulse can be deduced from the fact that, in the 
limit of a well-developed bow wave, all the electrons on a way of the laser pulse are pushed aside forming a multiflow 
motion owing to the collisionless nature of the low density plasma. The region void of electrons is positively charged 
and attracts the unpertubed electrons from aside. The transverse (radial) component of the electric field near the axis 
is proportional to the radius, E± — iT^n^er. It is easy to show that for the laser pulse waist smaller than i'Kca/ LOpe, 
where ujpe — 47rnee^/me)^/^ is the Langmuir frequency, the bulk electrons remain nonrelativistic. The electrons start 
to perform an oscillation with frequency f2 oc tOpe, determined by the geometry of the cavity. After one quarter of the 
oscillation period, 7t/2U, they reach the cavity axis. Therefore the characteristis shape of the cavity formed in the 
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FIG. 2: (color). 2D cross-sections of the electric field components (a) and Ey (d) and of the electron density (g) at 
t = 100 X 2ti/u. Vertical (b,e,h) and horizontal (c,f,i) ID cross sections of the E-field and density correspond to the positions 
shown by the dashed lines. In (a,d) the electron density curves are drawn by the black curves for rie/rieo = 0.75, 1.5. 



first wake wave period can be approximated by the curve 



r = i?cCOs(ri(i- x/vph)) 



(1) 



for — 7r/2 < Cl(t — x/vph) < 0, where Vph is the phase velocity of the wake wave equal to the group velocity, Vg, of the 
laser pulse. 



Comparing the kinetic energy of the expelled electrons. 



with their potential energy at the cavity periphery. 



eip — TTTieC^i?^, where Rc is the cavity radius, we find a condition of the bow wave formation: Rc < 2ca^^^ /ajp^ = 
{\ / '!T){ancr / neY^"^ ■ If we assume that the cavity transverse size is not greater than the laser pulse waist, Dias > Rc, 
then the condition of the bow wave formation can be rewritten in a stronger form as Dias < Aca^^'^ /ujpe. In the 
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FIG. 3: (color), (a) Magnetic field, B, around the cavities; laser frequency is filtered out. Red arrows for clockwise, blue - for 
counterclockwise vector B. On side panels ~ 2D cross-sections of the electron density curves for ne/rieo ~ 0.8, 1.4, 2, 2.6. (b,c) 
2D and ID cross section of By with electron density curves as in Fig. [2] Both frames for t — 100 x 2tv/uj. 
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FIG. 4: (color). Electron phase space projection onto (z,Pz) plane. 



opposite limit, Dias 3> 'ica^^'^ /ujpe, the momentum acquired by electrons from the laser pulse is not enough for 
escaping in transverse direction, so the bow wave is closely attached to the cavity boundary. 

As is well known, in one-dimensional (ID) case the laser pulse energy depletion length is given by Idep ~ hasi'^/'^pe)'^, 
where has is the laser pulse length, which comes from the balance of the laser energy and the energy left in the wake 
wave behind the pulse, when it is assumed that the electron energy is of the order of [T^. In the 3D 

configuration, when the laser pulse is tightly focused, Dias < 4:ca^^^ /wpe, so that the bow wave is generated, the 
kinetic energy of electrons being expelled in transverse direction is « amgC^, which results in the estimation for the 
depletion length: Idep ~ o-kasi^/^pe)'^ , i- e. by the factor a larger than in the ID case. 

In the LWFA concept [6 the energy of the electron accelerated by the wake field scales as 

£ « 2eA^7ph, (2) 

where Aip — (fmax ~ ^min is the potential difference in the wake field and 7p/i = (I — tip^j/c^)^^^^ is the Lorentz factor 
corresponding to the wake wave phase velocity. In the ID configuration the electromagnetic wave propagates with 
the group velocity, Vg, corresponding to 7g = (1 — = Iph — a^^^w/wpe in the limit a 3> 1. For the optimal 

conditions of the acceleration, eA(/3/meC^ — thus we obtain [T3] 

Eim^a^'^mec^ (3) 

^pe 
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Multi-dimensional effects modify the scaling law of the electron acceleration. In the cavity, formed behind a tighly 
focused laser pulse, close to the threshold of the bow wave formation, the electrostatic potential is mainly determined 
by the cavity transverse size, A(p « ~ irn^eDf^^/A:. For a narrow laser pulse its group velocity becomes strongly 
dependent on its waist, — Dias/^- The laser pulse waist, in its turn, depends on the pulse evolution, i. e. self- 
focusing and electron density cavity formation. We find that the energy of electrons accelerated in the wake wave 
generated by a sufficiently strong, a > eip/nieC^ — (7r/2)^(D;as/A)^(a;pe/Li;)^, and narrow, Dias ^ 4ca^/^/a;pe, laser 
pulse, does not depend on the laser pulse amplitude pTS] 

We note that here the condition for a implies E^^^^ < 87r~^a^(a;/a;pe)^rneC^- 

The shape of an outer boundary of the bow wave can be determined in the limit of a small density of electrons forming 
the bow wave, when these electrons can be treated as test particles. In this approximation, particles acquire initial 
momentum near the axis and then move in a potential of a positively charged semi-infinite wire of the cavity radius, 

i?c, propagating with velocity Vph- This motion is described by the Hamiltonian H — (rrilc^ + p^c 



2^2 I ^2^2 



1/2 



VphPx + ecr In ^+ {S,^ + y"^)^^'^ , for y > Re- On the particle trajectory in phase space the Hamiltonian is equal 

1/2 

to Ho = [mlc^ + p^qC^ + PyoC^) ~ VphPxO + ecr In [yo], where ^ = x — Vpfit, a w mieeR^ is the linear density of 
the wire, p^o and Pyo are the initial momentum components at y^ — Re- Neglecting the longitudinal momentum 
and assuming pyQ = arrive, we obtain ymax — Rc^xp [(1 -t- a^)^/^(2c/i?c'^pe)^] ■ Near the axis the outer boundary of 
the bow wave can be approximated by the straight line v±o(x — Vgt) + (vg — W||o)y = 0, where v\\q and wj^o £^re the 
components of the initial velocity of expelled electrons. The inclination of the line is determined by the expression 
tan^ — dy/dx = v±Q/{vg — wyo)- Extrapolating the approximation of the outer boundary of the bow wave by a line 
over the distance equal to the half-length of the cavity, we can estimate an addition to the cavity potential due 
to the bow-wave: tp^^'") « cr(l + ln[4 tan 6'/2i?J). 

Among various effects associated with the bow wave we notice the following. In the multi-species underdense plasma, 
on the ion time scale the bow wave structure becomes enriched by contribution from the ion dynamics. The bow wave 
has a complex structure, where modulations of lighter components develop faster, determining the coUisionless shock 
wave structure [3 ^] . The ion component of the bow wave is accelerated due to the laser ponderomotive potential 
and due to moving potential of the electron component of the bow wave. Transverse electric currents, associated 
with the bow wave from a laser pulse with strongly squeezed focal spot, should form large-scale quasi-static magnetic 
field and, due to instabilities, eventually should take a shape of separated jets. A portion of the transition radiation 
emitted at the plasma boundary interface [2(T can be attributed to a change of the electric charge density, associated 
with the bow wave propagation, especially in a non-uniform plasma. 

In conclusion, a tightly focused intense laser pulse propagating in underdense plasma, in addition to wake wave 
generation, excites bow waves and TWB waves, which spread the laser pulse energy in transverse direction. The 
effects of the bow and TWB waves excitation are essentially multidimensional. The bow wave increases the electric 
potential of the wake wave in comparison with the regimes where the bow wave can not detach from the cavity. The 
bow wave facilitates the transverse wave-breaking [TT] which causes the self-injection of electrons into the accelerating 
phase of the wake field. In its turn, the wave-breaking forms a transverse flow of electrons similar to the electron 
motion in the bow wave. 

Recent advances in LWFA [8j [12] reveal a tendency to use greater laser intensity and lower plasma density, aiming 
at higher energy of fast electrons. This makes the bow wave excitation inevitable. In such the regimes the multi- 
dimensional effects should modify the scaling laws of LWFA acceleration based on a one-dimensional approximation 
[H m [13] or on the assumption that the cavity size is determined by the laser pulse amplitude [^A^ . 
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